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We calculate the mass shift for the pion in a finite volume with renormalization 
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group (RG) methods in the framework of the quark-meson model. In particular, we 
investigate the importance of the quark effects on the pion mass. As in lattice gauge 
theory, the choice of quark boundary conditions has a noticeable effect on the pion 



\ mass shift in small volumes, in addition to the shift due to pion interactions. We 

in 

compare our results to chiral perturbation theory calculations and find differences 



due to the fact that chiral perturbation theory only considers pion effects in the finite 



Oh. volume. 



I. INTRODUCTION 

In the study of QCD, non-perturbative methods are essential in order to understand the 
connection between the high-momentum regime dominated by quarks and gluons, and the 
low-momentum regime described in terms of hadronic degrees of freedom. Lattice gauge 
theory is a method of great importance in this quest. Current simulations with dynamical 
fermions are limited to rather small lattice sizes and in some approaches to quark masses 
which are still large compared to the physical values. In addition to taking the continuum 
limit in which the lattice spacing is taken to zero, results from lattice calculations require 
extrapolation towards the chiral limit and the thermodynamic limit. Thus, in order to com- 
pare a result for an observable simulated in a small volume with the physical observable, it is 
essential to understand the finite volume effects. Apart from the application to lattice QCD, 
these finite volume effects are also interesting in their own right and worth investigating. 

The most important tool for extrapolations of lattice gauge theory results to small pion 
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masses and to large volumes is chiral perturbation theory (chPT) jl| 
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In particular for the chiral extrapolation to small pion masses 1. [l0|, llll Il2l LL3 



, and 



for the extrapolation to infinite volume for properties of the nucleon 3(, chiral perturbation 
theory describes the lattice results very well. 

In contrast to these applications, the finite volume shifts of the meson masses are less 
well described by chiral perturbation theory jbll I15I Q|. For the pion mass, the shifts 
predicted by chiral perturbation theory are consistently smaller than those observed in 
lattice simulations. We expect that chiral perturbation theory correctly describes finite 
volume effects for volumes that are sufficiently large so that the internal degrees of freedom 
such as quarks and gluons are unimportant 
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The discrepancies in current systematic 
llfil fl7 1 seem to indicate that there this 



investigations of finite volume effects 0, 
range has not yet been reached. 

One issue which cannot be addressed by chiral perturbation theory alone is the influence 
of the boundary conditions for the quark fields. While fermionic fields require anti-periodic 
boundary conditions in the Euclidean time direction, we are free to choose either periodic 
(p. b.c.) or anti-periodic (a. p. b.c.) boundary conditions in the spatial directions. In 
lattice calculations, this choice changes the finite size effects: In the investigation ^| of 
finite volume effects it was found that the choice of the boundary conditions for the quark 
fields has a direct influence on the size of the observed finite volume shifts jlfll . and 
an explanation in terms of quark effects was proposed, for both quenched and unquenched 
calculations. Such effects cannot be captured by a description in terms of pion fields only. 
It has been shown by Gasser and Leutwyler that the low-energy constants in the chiral 
perturbation theory Lagrangian remain unchanged from their values in infinite volume if 
one considers QCD in a finite Euclidean volume, provided the same anti-periodic boundary 
conditions as in the temporal direction are chosen as well in the spatial directions for the 
quark fields jjj. This leaves open the possibility that a change of boundary conditions for 
the quark fields might in fact lead to a change in the finite volume behavior. 

To motivate further our interest in the influence of boundary conditions, in Fig. we 
present an example of a lattice calculation in the quenched approximation from the ZeRo- 
collaboration [lfij . Shown is the shift of the pion mass m^^L) in finite volume relative to 
the value in infinite volume m n (oo) as a function of 777^(00) ■ L where L is the volume size. 



Surprisingly, these results show a dropping pion mass for intermediate volume sizes in a 



3 



R m (L) = (m^-mj/rr^ 




k = 0.1340 
0.1345 
0.1350 



J i- 



9 10 11 12 



FIG. 1: The pion mass shift R[m n (L)] = (m n (L) — m ^jo o)) I m^jco) as a function of m n (oo) ■ L, 
obtained in a quenched lattice calculation, from ref. |l5|. Shown are results for three different 
values of the quark mass, determined by k. The solid lines show the corresponding predictions 
from chiral perturbation theory. 
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FIG. 2: The pion mass m 7r (L) as a function of box size L, obtained in a lattice simulation 
with two flavors of dynamical Wilson fermions, from ref. Results for pion masses = 

643 Me"VVi90 MeV, 419 MeV (circles, squares, diamonds) are compared to results from Ltischer's 
formula ^| with input from chPT up to NNLO order !4|. For details, see ref. |l(J]. For the smallest 
pion mass, a drop similar to the one in Fig. ^ can be observed. 
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region where the standard chiral perturbation theory result (indicated by the solid lines in 
the figure) predicts only a very weak volume dependence. This behavior would be unexpected 
from pion effects alone. In addition, finite volume effects from chiral perturbation theory 
are predicated on the presence of a "pion cloud", which in turn requires the presence of 
sea quarks Like chPT, the present work is not directly concerned with the quenched 



approximation, which requires its own low-energy effective theory 21, 22]. Although the 
quenched calculation in Fig. ^ shows the pion mass drop in a very distinct fashion, similar 
effects for the meson masses are also seen in studies of finite size effects with dynamical 
quarks Q] . In Fig. |21 we show results for the pion mass as a function of the volume size 
from a lattice calculation with two dynamical flavors of Wilson fermions r| . Results are 
given only for few volume sizes, but they also show a decrease of the pseudoscalar mass in 
small volumes. 

In this paper, we investigate the volume dependence of the pion mass in the quark meson 
model using renormalization group (RG) methods. In particular, our purpose is to extend 
our previous work [23| and to investigate the influence of different boundary conditions for 
the fermionic fields on the finite volume effects for low-energy observables such as the pion 
mass and the pion decay constant. Since our model includes dynamical meson fields, and a 
dynamical breaking of the chiral flavor symmetry according to SU(Nf) x SU (Nf) — > SU(Nf) 
for Nf = 2 flavors of quarks, our results are applicable to unquenched lattice calculations 
with two dynamical quark flavors. 

The quark-meson model cannot predict the volume dependence of pion mass and pion 
decay constant exactly. It is not a gauge theory and thus has neither gluons nor quark 
confinement. At moderate energies, below the hadronic mass scale, unconfined constituent 
quarks appear instead of baryonic degrees of freedom. On the other hand, the model ha s 
been rather successfully employed in the description of the chiral phase transition (24 , 
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the linear sigma-model with quarks are compatible with those of 



The low-energy couplings o 
chiral perturbation theory 2jj. As we have previously shown 23], for small pion masses and 
large volumes our results for the volume dependence agree with those of chiral perturbation 

nn 

theory |4J, |£J , if we apply anti-periodic boundary conditions for the quark fields in the spatial 
directions. In this paper, we investigate the effect of different boundary conditions for the 
quark fields on low-energy observables, namely the pion mass and pion decay constant, in 
more detail. 
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In spite of the shortcomings of our model, we believe that the current approach can shed 



light on lattice results regarding the volume dependence of the pion mass 
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ljj. While 



the actual mechanism in QCD may be different due to the presence of color interactions, the 
approach employed in the present paper gives a possible explanation for the apparent drop 
in the pseudoscalar pion mass in small volumes observed in 14, 3,0], which precedes the 



rise of this mass due to chiral symmetry restoration in extremely small volumes. 

The paper is organized as follows: In section |H1 we briefly introduce the model and the 
renormalization group equations which govern the RG flow in a finite volume. In section ITTT1 
we solve the flow equations numerically and in section IIVI we present the results for the 
volume dependence of the pion mass. A comparison to lattice results and our conclusions 
are found in section El 



II. RG-FLOW EQUATIONS FOR THE QUARK-MESON MODEL 

As motivated in the introduction, we will use an 0(4)-invariant linear a- model with 
Nj = 4 mesonic degrees of freedom (cr, tt), coupled to Nf = 2 flavors of constituent quarks 
in an SU(2)l ® SU(2)f> invariant way. This model does not contain gluonic degrees of 
freedom, and it is not confining, but it is an effective low-energy model for dynamical 
spontaneous chiral symmetry breaking at intermediate scales of k < hjjy. The ultraviolet 
scale Kjjv ~ 1-5 GeV is determined by the validity of a hadronic representation of QCD. At 
the UV scale Auv, the quark-meson-model is defined by the bare effective action 

T Auvl0} = J d 4 x {q($ + gm c )q + gq(cr + ir ■ 7rj 5 )q + l(d^) 2 + U Auv ((p)} (1) 

with a current quark mass term gm c which explicitly breaks the chiral symmetry. The 
mesonic potential is characterized by two couplings: 

UK v A4>) = \rn 2 uv 4> 2 + \\uv{4> 2 ) 2 . (2) 
In a Gaussian approximation, we obtain the one-loop effective action for the scalar fields 0, 
m = r Avv [<P] - Tr log (if [0]) + ^Tr log (if [</>]) (3) 

where if [0] and if [0] are the inverse two-point functions for the bosonic and fermionic 
fields, evaluated at the vacuum expectation value of the mesonic field 0. We consider the 
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effective action T in a local potential approximation (LPA), where the expectation value is 
taken to be constant over the entire volume. In order to regularize the functional traces, 
we use the Schwinger proper time representation of the logarithms. A scale dependence is 
introduced through an infrared cutoff function 



t|/.(r*«) 



irk') 



2\a+l e -rk 2 



(4) 



T(a + 1) 

which regularizes the Schwinger proper time integral. By replacing the bare coupling in 
the inverse two-point functions with the scale-dependent running couplings, we obtain a 
renormalization group flow equation for the effective potential in infinite volume for zero 
temperature: 



d 

k^rU k (a, vf 2 , T — > oo, L — > oo) 
ok 



fc2(a+l) 



AN c N f 



16a(a - 1)tt 2 
1 



+ 



(k 2 + M 2 (a, tt 2 ))"- 1 
Nj-1 



(5) 



(A; 2 + M2(a,7T 2 ) a - 1 (k 2 + M 2 (a, tt 2 ))"- 1 

In infinite volume, a = 2 is the lowest possible integer value we can choose in order to be able 
to perform the Schwinger-proper time integration. Note that in LPA, the effective action 
reduces to the effective potential through the relation 



T fe [0] = J d 4 xU k (a,n 2 ). 



(6) 



Due to the fact that we allow for explicit symmetry breaking, the effective potential 
becomes a function of o and tt 2 . The meson masses are the eigenvalues of the second 
derivative matrix of the mesonic potential: 



2U^+A tt 2 CW + U aa + ^(2 + 4 n 2 U^2 - U aa ) 2 + 16 vf 2 U 2 



1 



M 2 



2^2+4: vf 2 U^2 + U aa - y/(2 U^2 + 4 FCW - U aa ) 2 + 16 vf 2 U 2 



(7) 



For 7? 2 = 0, the masses of the three pion modes are degenerate. The symmetry breaking 
terms in the a-direction do not affect the 0(3)-symmetry of the pion subspace, so that the 
pion fields appear only in the combination n 2 in the eigenvalues. The constituent quark 
mass is given by 



Ml = g 2 [(a + m c f + tt 2 ] 



(8) 
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To derive renormalization group flow equations in a finite four-dimensional volume L 3 x T, 
we replace the integrals over the momenta in the evaluation of the trace Q by a sum 

dft-..- T E ••• ( 9 ) 

rii=— 00 

We are free in the choice of boundary conditions for the bosons and fermions in the space 
directions. However, in the time direction, the boundary conditions are fixed by the statistics 
of the fields. Adopting the language of lattice literature, we use T to denote the length of 
the finite volume box in Euclidean time direction . Then the Matsubara frequencies take 
the values 

2 ™o , (2n + l)7r 
^n = ~y~ and u no = . (10) 

for bosons and for fermions, respectively. In the following we use the short-hand notation 

, 4vr 9 :! ' ' >; 



P P 



4f 2 A , , 2 4vr 2 ^ / IV , . 

JTl^ n i aIld Pap = JY 1^ K + 2 J (U) 

i=l i=l ^ ' 



for the three-momenta in the case of periodic (p) and anti-periodic (ap) boundary conditions. 
The flow equation corresponding to eq. (jHJ) is 

d „ , ^ 2 P( a+1 ) >r^>r^ / AN c N f 



n n \ n 1 1 ""<2 



1 \ 

U (^ 2 + <+^ + ^(a,7r2))^J- (12) 

The sums in eq. f!12|) run from — oo to +00, where the vector n denotes (n 1; n 2 , 72,3). For 
both finite and infinite volume, we choose a = 2 for the cutoff function. For a volume with 
infinite extent in time direction T — ► 00, we perform the sum over the Matsubara frequencies 
analytically 28]: 

h d_ n( . 2t n-A*!.W *Ml 

iV?=4 



2 (p +p 2 + M 2 (a;7f2))5 /2j ( 13 ) 



Note that we employ both flow equations (J12|) and (|13|) for our numerical calculations in 
the next section. We would like to make one comment: The insertion of the regulator 
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function in the Schwinger proper-time integral is not necessary to regularize the infrared 
regime, since finite volume calculations are already infrared finite. However, if we keep the 
volume fixed, the insertion of the regulator function is needed to integrate out the quantum 
fluctuations in a controlled way. 

In order to solve the partial differential equations (jHJ and ((T2j) , we project these flow 
equations on the following ansatz for the mesonic potential 2^ |: 

f4(a,7f 2 ) = E ^(k)^-Mk)Y(cr 2 + n 2 -a (k) 2 y. (14) 

i=0 j=0 

Performing such a projection, we get, in principle, an infinite set of coupled first-order 
differential equations. To solve this set of equations, we have to truncate the ansatz at some 
power in (er — cr (k)) and (cr 2 + n 2 — a (k) 2 ). In this paper, we use N a = 2. The resulting 
finite set of flow equations can be solved straightforwardly in a numerical calculation. 

III. CALCULATION 

We have solved the RG flow equations numerically and will present the results for the 
volume dependence of the pion mass and the pion decay constant in the following section. 
First we discuss some details about the numerical evaluation and the determination of the 
coefficients of the ansatz for the potential eq. ()14j) at the UV scale. At the ultraviolet 
cutoff scale Ajjv, the meson potential can be characterized by the values of the couplings 
muv and Xuv i n eq- 0- All other coefficients in the ansatz eq. ()14)) are set to zero. In 
order to solve the flow equations for the effective potential, we truncate the ansatz for the 
potential as discussed in the last section. In the present work, we expand the potential up 
to mass dimension four in the mesonic fields. Furthermore, we have to specify a value for 
the current quark mass m c , which controls the degree of explicit s ymm etr y b reaking. The 
Yukawa coupling g does not evolve in the present approximation 24, 0, 13] • We choose 
g = 3.26, which leads to a reasonable constituent quark mass of M q = gif-K + m^ ~ 310 MeV 
for physical values for the pion decay constant f v = 93 MeV and the current quark mass 
gm c = 7 MeV. 

In table HI we summarize the three parameter sets which we use in obtaining our results 



for pion masses of 100, 200 and 300 MeV, see also ref. (21 



. In our comparison of different 



boundary conditions, we use the same parameter sets to obtain results for either periodic 
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A uv [MeV] 


m uv [MeV] 




gm c [MeV] 


U [MeV] 


[MeV] 


1500 


779.0 


60 


2.10 


90.38 


100.08 


1500 


747.7 


60 


9.85 


96.91 


200.1 


1500 


698.0 


60 


25.70 


105.30 


300.2 



TABLE I: Values for the parameters at the [TV-scale used in the numerical evaluation. We 
determine these parameters by fitting in infinite volume to a particular pion mass m 7r (oo) and the 
corresponding value of the pion decay constant f n (po), taken from chiral perturbation theory. In 
our notation, the physical current quark mass corresponds to gm c . We use g = 3.26. 

or anti-periodic boundary conditions for the fermions. We determine these UV parameters 
by fitting to the pion mass m n (oo) and to the corresponding pion decay constant /^(oo), 
which is taken in infinite volume from chiral perturbation theory ^J. We then evolve the 
RG equations with these parameters to predict the volume dependence of fn(L) and m^L). 
From table 12 we can read off that the pion mass is primarily determined by the value of the 
current quark mass, which controls the explicit symmetry breaking. To obtain the correct 
value for the pion decay constant for a given pion mass, the meson mass at the UV scale 
muv has to be decreased from 780 MeV to 700 MeV when the pion mass increases from 
100 to 300 MeV. We use the same value for the four-meson-coupling Xuv f° r all values of 
the pion mass and pion decay constant considered here. The possible values of the current 
quark mass are limited by the requirement that all masses, in particular the sigma-mass, 
must remain substantially smaller than the ultraviolet cutoff Auv ~ 1500 MeV of the model. 
We have checked that our results are to a large degree independent of the particular choice 
of UV parameters: different sets of starting parameters give the same volume dependence, 
provided that they lead to the same values of the pion mass and pion decay constant in 
infinite volume. 

We use the result of chiral perturbation theory for the dependence of the pion decay 
constant on the pion mass to facilitate the comparison between the quark meson model and 
chPT. However, it is possible to get the correct behavior of the pion decay constant as a 
function of a single symmetry breaking parameter with renormalization group methods, as 
was shown in infinite volume 
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For anti-periodic boundary conditions for the fermions, we have previously investigated 

n 

the dependence of our results on the cutoff scale Ajjv .23] and found it to be small for 
light pions, and only moderate for the largest pion mass we considered here. We argued 
that such a trend towards a stronger cutoff dependence for larger pion masses was to be 
expected, since the existence of a cutoff becomes more relevant for heavier mesons. This 
analysis still pertains to the results with anti-periodic boundary conditions presented here. 
For a choice of periodic boundary conditions, however, the cutoff dependence of the results 
is somewhat more pronounced, in particular for small volumes when the Euclidean time 
extent is kept large. Varying the cutoff between 1.5 GeV and 1.1 GeV for a pion mass of 
77^(00) = 300 MeV, we find that the largest variations are of the order of 5 — 6% of the 
pion mass, and take place in a volume range of L = 0.5 — 1.0 fm, depending on the exact 
ratio T/L of time and space extent. As we argue below, this is mainly due to effects on 
the quark condensation: for periodic boundary conditions, a larger UV cutoff allows for the 
build-up of a larger condensate in finite volume, since for any given volume, a larger number 
of momentum modes 2n\n\/L remain below the cutoff and contribute. In a volume region 
where the quarks dominate the finite volume effects, a certain cutoff dependence of these 
effects is therefore expected in this model. 

The sums over the momentum modes in the flow equations cannot be performed analyt- 
ically, therefore we have to truncate the sums at a maximal mode number N max = max|n|. 
With this truncation, we introduce an additional UV cutoff in our calculation. In order to 
guarantee that this cutoff does not affect our results we require 

-j-N max > A uv . (15) 

We have to take care that this relation is well satisfied since we are using a "soft" cutoff 



function. We have checked the dependence of the results on N max in [23j and found that it 
is sufficient to use N max = 40 for Ajjv = 1.5 GeV and volumes up to L = 5 fm, which we 
will also use for the calculations in this paper. The numerical evaluation of the sums over 
the momentum modes simplifies significantly if we take the box sides as integer multiples of 
some length scale Lq, such that L = ulLq and T = titLq. Therefore we restrict ourselves to 
this case. Below, we show that the results for the low-energy observables strongly depend 
on the ratio T/L = n^/n^. 
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IV. RESULTS 

We have calculated the pion mass shift 

R[m n {L)\ = — (16) 

with both choices for the fermionic boundary conditions for three different pion masses, 
m n (oo) = 100, 200 and 300 MeV, and for infinite (T/L — > oo) as well as for finite extent of 
the Euclidean time axis with different ratios T/L = 3/1, 3/2, 1/1. 

In Fig. H3 we show the results for the pion mass shift with periodic boundary conditions 
as a function of the box size L. The three panels show the results for the three different pion 
masses we investigated, and the curves are labeled with the ratios T/L. The main new and 
surprising observation is that in this case, for certain volume ranges, the mass of the pion in 
the finite volume can be lower than in infinite volume. In particular, this is the case for pion 
masses m n (oo) > 200 MeV, ratios T/L > 3/2, and volumes smaller than 2fm: R[m n (L)] 
takes on negative values and develops a minimum. This can be seen in the lower two panels 
of Fig. 01 Secondly, we note that this minimum in the mass shift becomes deeper for larger 
pion masses m n (oo), and the corresponding larger values of f n (oo). For m n (oo) = 300 MeV, 
the pion mass shift reaches down to approximately i2[m„.] = —0.14 at L = 0.7 fm. 

In Fig. HJ we compare the results for the pion mass shift with periodic (p. b.c.) and anti- 
periodic (a. p. b.c.) boundary conditions for the fermion fields, for the ratios T/L = 3/2 and 
T/L = 1/1. Clearly, employing p. b.c. lowers the relative mass shift R[m n (L)], compared 
to using a. p. b.c. The differences become larger in smaller volumes, for larger pion masses 
m % (oo), and with increasing ratios T/L. As we have seen, if the length of the box in the 
Euclidean time direction is taken to infinity, for large pion masses the pion mass can be 
smaller in finite than in infinite volume, so that the finite volume shift becomes negative. 

Although at first a surprising result, this shift to smaller pion masses can actually be 
explained in the framework of the quark-meson model and its mechanism of chiral symmetry 
breaking. In order to show this, we resort to a version of the model that is simplified 
compared to our ansatz (|14p. but still contains the same essential structure. In this model, 
for a fixed symmetry breaking parameter gm c , the pion mass is completely specified by the 
scale-dependent order parameter a (k,L), and by the values given at the UV scale for the 
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FIG. 3: Results for the pion mass shift i?[m„-(L)] = (m 7r (L)— m 7r (oo))/m vr (oo), in a finite Euclidean 
volume of size V = L 3 x T, for periodic boundary conditions. The ratio of T/L for the different 
curves is given in the figure. We show the results for pion masses of m^-(oo) = 100,200,300 MeV 
(identified in the figure). 
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coupling g and the meson mass mfjy. According to 0, |^1|, it is 

M^,L) = ^f%. (17) 
cr {k,L) 

For periodic boundary conditions, the "squeezing" of the quark fields in a small finite volume 
leads to an increase in the chiral quark condensate, before a further decrease of the volume 
size induces a restoration of chiral symmetry. Following eq. (|17jl. the increase in the order 
parameter leads in turn to the observed decrease in the pion mass. 

The intermediate increase in the order parameter with the decreasing volume size can be 
explained more rigorously from the flow equations. Since this increase occurs in volumes 
that are already quite small, the flow is dominated by the zero-momentum modes and it is 
sufficient to analyze the contributions of these modes. 

The zero mode contribution to the flow equation from quarks and mesons is for purely 
periodic boundary conditions in spatial directions given by 



d 



k 2( - a+ V ( AN c N f 



o 



TL 3 V (k 2 + ui + Mg(k,L,a,Tf 2 ) 2 Y a +^ 
Nj-1 1 \ 

~(k 2 + M n (k,L,a, vr 2 ) 2 )^ 1 ) ~ [k 2 + Mjk, L,a, vr 2 ) 2 )^ 1 ) ) (18) 

where Uq = (±7r/T) 2 corresponds to the value of the two Matsubara frequencies closest to 
zero. The prefactor 1/L 3 diverges for L — > for all momentum modes, but enhances only 
the zero modes: For the non-zero momentum modes, the enhancement is canceled and they 
are in fact strongly suppressed, which is due to the factors 1/L 2 of the momentum terms in 
the denominators. If we scale T proportional to L, because of the Matsubara frequencies 
this suppression occurs also for the lowest fermionic terms, although it is much weaker. The 
result of this competition between suppression and enhancement for the fermions depends 
on the ratio T/L. 

We first consider exclusively the contributions of the fermionic zero modes, which exist 
only for periodic boundary conditions: 

d n 1 F k 2( - a+ ^ AN Nf 

k m U k (c,t\T,L)^ =-_ .2- WT1¥T ^J-^ (19) 

This truncation to the fermionic contributions only is equivalent to the leading term of a large 
iV c -approxiination, as it is shown in j^. In principle, eq. (|T9|) can be integrated analytically, 
since the constituent quark mass, given by M 2 (a, 7? 2 ) = g 2 [(cr + m c ) 2 + 7? 2 ], does not depend 
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on any scale-dependent quantities. The result shows that the zero mode contributions to the 
potential as a function of the expectation value are repulsive for small values. Consequently, 
these contributions increase the expectation value ao(k,L) and thus the value of the pion 
decay constant. Since these zero- momentum contributions are enhanced for small volumes, 
this explains the increase in the expectation value. 

Alternatively, this can be understood in more detail by a direct analysis of the zero mode 
contributions to the flow equation for the minimum <Jo(k, L) of the potential. Since the flow 
equation for ao(k, L) is obtained from the minimum condition 

= *o(k, L),tt 2 = 0, L, T) = 0, (20) 

it is determined by the flow of the potential. As we have seen in our analysis above, the 
fermionic contributions tend to increase the absolute value of the minimum ao(k,L), while 
the mesonic contributions tend to decrease it. Thus, we can perform this analysis entirely 
by considering the zero mode part of the potential flow given in eq. (|18|). 

The renormalization scale k controls the momenta of the quantum fluctuations that are 
integrated out. As soon as this momentum scale drops below the mass of one of the degrees 
of freedom, that particular field can no longer contribute to the RG evolution of the running 
couplings: it decouples from the RG flow. We restrict the discussion here to scales k < m a , 
where the sigma meson has already decoupled. 

With periodic boundary conditions, the finite box length in the Euclidean time direction 
T is the only scale which affects the zero modes. The scale n/T is in competition with the 
renormalization scale k, and if k drops below this scale, the lowest Matsubara frequency 
uq = tt/T acts as a cutoff and stops that part of the evolution which is driven by the 
quark fields. If T is sufficiently small, this happens already above the scale at which chiral 
symmetry breaking sets in. In that case, condensation of the quark fields is prevented, and 
the constituent quark mass remains small. This means that m n (k — > 0,L) remains large 
and that i?[m^(L)] is large and positive. This is illustrated by the results for T/L = 1/1 
and small L in Fig. El 

The situation is different for large values of T/L > 3/2. Here, the additional scale set 
by 1/T plays a less important role and becomes relevant only for much smaller volumes. 
In this case, quarks build up a large condensate. According to eq. ([17)1 . this increase in 
the chiral condensate leads to a decrease of the pion mass, which is visible in Fig. El for 



15 



T/L > 3/2, m 7r (oo) > 200 MeV, and L > 0.8 fm. For large values of T/L, the decrease in 
the condensate for small volumes cannot be explained by the presence of the cutoff n /T for 
the quark fields alone. There is an additional mechanism that decreases o"o in such a way 
that chiral symmetry is broken less strongly. For very small volumes, the pion contributions 
in eq. (jTHj) dominate the flow of ao- Even for a large ratio T/L, this leads to a decrease in 
a and the observed rise in W^m^L)] for small L. 

For anti-periodic boundary conditions, we do not find any decrease of R[m n (L)} with 
decreasing finite volume size L for any value of T/L, as can be seen in the comparison 
in Fig. HJ As we have argued in 2^, in our RG approach with anti-periodic boundary 
conditions, two effects are responsible for the finite volume behavior: effects due to the quark 
condensation, and effects due to light pions which appear after the chiral condensate has 
been built up by the quark fields. In contrast to the case of periodic boundary conditions, 
for anti-periodic boundary conditions the formation of the quark condensate is strongly 
suppressed by the lowest possible momentum for the fermions, which is y/3ir/L, see eq. ([lip, 
and acts as an infrared cutoff. Consequently, for small L, fewer modes contribute to the 
chiral condensate. If in addition T/L is small, the condensate decreases further and we 
observe a larger mass shift R[m n (L)}. 

Finally, in Fig. |5] we compare our results for the pion mass shift to the results of chi- 
ral perturbation theory. (Note that Fig. |5] has a logarithmic scale, whereas Figs. E2 and 0] 
have linear scales.) We present results for different pion masses from RG calculations with 
both periodic and anti-periodic boundary conditions for the fermions, and from chiral per- 
turbation theory 0, E|. For the chPT results, the pion mass shift is calculated with the 
help of Liischer's formula 18], which relates the leading corrections of the pion mass in 
finite Euclidean volume to the 7T7r-scattering amplitude in infinite volume. The sub-leading 



corrections drop as 0(e~ mL ) with fa > a/3/2 m n . Using a calculation of the 7T7r-scattering 
amplitude in chPT to three loops (nnlo) as input for Liischer's formula, the authors of ref. $ 
obtain a correction above the leading order, which is then added to the one-loop result of 
Gasser and Leutwyler j^j]. Liischer's original approach only considers the periodicity of 
pion propagators in finite volume as an invariance under a shift by L. More recently, in 
ref. p| this has been improved to account for the fact that these propagators are actually 
invariant under shifts by nL with arbitrary ft. The result is a Liischer formula resummed 
over n, which is very similar to the original one. The finite volume shift for the pion mass 
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FIG. 4: Comparison of the pion mass shift in finite volume R[m w (L)] = (m 7r (L) — m 7r (oo))/m 7r (oo) 
for the two choices of fermionic boundary conditions. Open symbols denote results for anti-periodic, 
solid symbols for periodic boundary condition. The size of the volume is V = L 3 x T, the ratios 
of T/L for the different curves are given in the figures. We show results for pion masses of 
m^oo) = 100, 200, 300 MeV (identified in the figure). 
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is significantly increased by this resummation. In 



we have carefully compared our RG 



results with anti-periodic boundary conditions to the chPT results from |4j. Here, we use the 
improved results from (| for the comparison. The RG results are still consistently above the 
results from chPT. Liischer's approach becomes an increasingly better approximation with 
increasing pion mass for a given volume size. The decreasing differences between the chPT 
results and the RG results with increasing pion mass are compatible with this estimate. For 
large volumes, the mass shift is completely controlled by pion effects and drops as e _m " rL , 
so that both the RG and the chPT results have the same slope in the logarithmic plot. For 
the entire volume range shown in Fig. the RG and chPT results apparently differ only 
by a factor which is almost independent of the volume size. For 77^(00) = 300 MeV, the 
chPT and RG results agree within errors. For small volumes, however, the RG approach 
has the advantage that it can be extended to describe the transition into a regime with 
approximately restored chiral symmetry, where the chiral expansion becomes unreliable. 

The mesonic degrees of freedom are less affected by the ratio T/L. The upper curve in 
Fig. |3] represents RG calculations with anti-periodic boundary conditions and T/L = 1/1, 
which gives a larger R[m n (L)\ compared to the lower curve corresponding to T/L = 00. 
Fluctuations due to the light pions yield a decrease of the condensate and explain the increase 
of -Rfm^L)] for larger volumes. In particular for small pion masses (m^ = 100 MeV) and 
the ratio T/L = 1/1, the results with periodic and with anti-periodic boundary conditions 
overlap over a wide volume range. From our analysis, for sufficiently small values of 77^(00) 
this is expected in the region where pion dynamics dominate. Because of this, the slopes of 
the curves are very similar. The deviations between results at the same, fixed ratio T/L that 
differ only in the choice of boundary conditions become larger for increasing pion masses 
772^(00) and decay constants /^(oo). This indicates that fermionic effects are increasingly 
important. Evidence for this is also the observation that the results for the pion mass shift 
with periodic boundary conditions have a smaller slope, compared to the results with anti- 
periodic boundary conditions, and also compared to those of chPT. The reason is that the 
cutoff scales are different: for periodic boundary conditions, the lowest fermion momentum 
mode is given by the lowest Matsubara frequency uq = n/T, and not determined by V3n/L 
as for anti-periodic boundary conditions. In particular for large values of T/L, this explains 
that the finite volume mass shift will be much larger for anti-periodic boundary conditions. 
For small volumes, we thus find the importance of quark effects confirmed by the dependence 
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on the boundary conditions. But since pion effects dominate for larger volumes, the results 
of chPT and of our RG approach converge for large L. 

V. COMPARISON TO LATTICE RESULTS AND CONCLUSIONS 

Apart from the general interest of finite volume effects, the main motivation for our 
current investigation is its possible application to lattice gauge theory. At present, most 
lattice calculations are performed in volumes of the order of L = 2 — 3 fm. In recent 
systematic studies of finite volume effects done with Wilson fermions, the lightest pion 
masses are of the order of m n = 400 - 500 MeV Q, Q HQ, H, Q- 

With staggered 

fermions, pion masses as low as 250 MeV have been realized j^Sj . Simulations with fermions 
with good chiral properties such as domain-wall or overlap fermions have been done with 



pion masses as low as 180 MeV in the quenched approximation [36J, |37( and as low as 
360 MeV with two fully dynamical flavors Because the finite volume effects depend 

on the mass of the lightest field, they become more severe for smaller pion masses. Thus, 
the better the statistical accuracy of these calculations, the more important it becomes to 
understand finite size effects and to control the finite size extrapolation. 

Our model incorporates chiral symmetry and can still be used in the vicinity of the point 
where chiral symmetry is restored. Finite volume effects should therefore be captured as 
far as they relate to chiral symmetry breaking. But it is not a gauge theory, there are 
no gluons, and consequently the constituent quarks in this model are not confined. There 
is no guarantee that the same mechanisms apply as in QCD. Since the model contains 
dynamical meson fields and chiral symmetry is broken in the usual way, our results can 
only be compared directly to those of unquenched lattice calculations with two dynamical 
quark flavors, where normal chPT is also applicable. However, qualitatively our arguments 
regarding the quark condensate may also have implications for quenched simulations, since 
a similar mechanism may apply. 

Extrapolations to infinite volume using chiral perturbation theory are extremely success- 
ful in the description of the volume dependence of nucleon properties, such as for example 
the nucleon mass 

a a. 

However, as far as meson masses are concerned, the finite vol- 
ume mass shifts observed on the lattice deviate from the predictions of chiral perturbation 
theory. This holds also Q| for Liischer's approach ^], which only takes pion effects into 
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account as well. Generally, the predicted mass shifts are much smaller than the observed 
ones 14,15, r| 19 ] . The inclusion of higher orders in the chiral expansion ^| and a summa- 
tion of additional contributions in Luscher's expression (f| increase the size of the predicted 
mass shifts and decrease the distance to the RG results. For the physical values of the pion 
mass and the pion decay constant, chiral perturbation theory can be applied for volume 
sizes L ^> lfm, but a priori it is impossible to say how large exactly the volume has to be, 
according to ref. Ultimately, this question can only be answered by lattice calculations. 

Comparing our results to those from chiral perturbation theory, we find agreement for 
larger pion masses, provided we impose anti-periodic boundary conditions on the fermionic 
fields. As expected, the differences increase for very small volumes, where chiral symmetry 
restoration becomes important. For periodic boundary conditions, large pion masses, and a 
large ratio T/L, the mass shifts in small volumes behave differently from those of chPT. In 
particular, even in volume sizes as large as L = 2.5 fm and for T — > oo, the results for the 
relative shift of the pion mass can as much as double under a change of boundary conditions, 
for example from R[m n (L)] = 0.0488 with periodic to R[m n (L)] = 0.1037 with anti-periodic 
boundary conditions for a pion mass of m n = 100 MeV. In Table |H] we give bounds on 
the minimum size of the volume that are necessary to keep the finite volume mass shift 
smaller than 10% resp. 1%, calculated for periodic and antiperiodic boundary conditions 
in the RG approach, and for comparison from the NNLO chPT calculation of ref. In 
general, periodic boundary conditions allow to achieve the same accuracy with regard to 
finite volume effects with smaller volume sizes than anti-periodic boundary conditions. 

For periodic boundary conditions, our results reproduce the qualitative behavior of the 
lattice results, but clearly differ from chPT. For anti-periodic boundary conditions, they 
largely agree with chPT. This suggests that for anti-periodic boundary conditions, the ef- 
fective low-energy constants relevant for the finite volume effects of our observables agree 
with those of chPT, in agreement with the argument by Gasser and Leutwyler for QCD 0. 
However, the differences then might imply that the low-energy constants change for periodic 
boundary conditions. 

The issue of finite volume effects has been addressed in several lattice studies llj, Il5l 
The pion mass shift R[m n (L)\ calculated by the ZeRo collaboration , which 
is shown in Fig. Q actually becomes negative and has a minimum at small volume sizes. 
Although this negative shift is small, the result seems to be significant. The minimum is 
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m„-(oo) 



R[m n (L)] < 0.1 



R[m^(L)\ < 0.01 



100 MeV 



RG, ap. 
RG, p. 
chPT 



L > 2.523 fm 
L > 1.351 fm 
L > 2.187 fm 



L > 4.381 fm 
L > 4.259 fm 
L > 3.785 fm 



200 MeV 



RG, ap. 
RG, p. 
chPT 



L > 1.736 fm 
L > 0.5 fm* 
L > 1.639 fm 



L > 2.842 fm 
L > 1.888 fm 
L > 2.653 fm 



300 MeV 



RG, ap. 
RG, p. 
chPT 



L > 1.359 fm 
L > 1.022 fm 
L > 1.339 fm 



L > 2.213 fm 
L > 1.911 fm 
L > 2.104 fm 



TABLE II: Bounds on the minimum size of the volume V = L 3 x T for T — > oo such that the finite 
volume pion mass shift R[m n (L)] is < 0.1 or < 0.01, for different values of the pion mass 771^(00). 
RG results are given for anti-periodic and for periodic boundary conditions, chPT results are those 
in NNLO obtained in |6j. Note that for periodic boundary conditions and for m n = 200 MeV 
and m n = 300 MeV, the bounds are set by a decrease of the pion mass. *For = 200 MeV 
and periodic boundary conditions, the bound R[m n (L)] < 0.1 is satisfied in the full volume range 
described by our model (cf. also Fig. |HJ). 

most pronounced for small quark masses (at k = 0.1350). The position of the minimum 
corresponds to m^L = 3.5 or L = 1.264 fm with T/L = 2.25. The results were obtained in 
the quenched approximation with periodic boundary conditions for the quark fields. Similar 
observations have also been made in 14 , ^| , where the simulations were performed with 
dynamical Wilson quarks. 

In our calculation, such a decrease in the pion mass is reproduced if we choose periodic 
boundary conditions for the quarks. The minimum appears for large pion mass m^oo) = 
300 MeV, T/L > 3/2 and L = lfm, cf. Fig. |H1 Our model suggests a mechanism for 
the appearance of this minimum, which may be the same mechanism as on the lattice. 
In contrast to our findings, however, the decrease of the pion mass in finite volume seems 
to be larger for smaller infinite-volume pion mass. For lattice calculations, several other 
mechanisms for finite volume mass shifts have been suggested, from an interaction of hadrons 
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with their mirror states on a periodic lattice |39j to effects on quark propagation related to 
a breaking of the center symmetry of the gauge group |ig| | . 

The influence of boundary conditions for sea and valence quarks in lattice simulations 
was also studied by Aoki et al. [19]. They find that periodic boundary conditions lead 
to a lower mass shift than anti-periodic boundary conditions (see table III of [19]). This 
finding is in agreement with our results, as can be seen in Figs. 0] and 03 The actual pion 
mass on the lattice is very high (> 1 GeV). Different choices for the boundary conditions of 
sea and valence quarks make it possible for the authors to establish a connection between 
the mass shift and the expectation value that Polyakov loops acquire in the presence of 
sea quarks. They relate the large increase of the pion mass observed for small lattice size 
to the restoration of chiral symmetry. This is illustrated by their results for the chiral 
condensate (Fig. 10 of which decreases strongly in small volumes. In the same figure, 
the condensate may increase for intermediate volume size, which would be similar to the 
behavior of the order parameter seen in our simple model. We agree that the mass shift 
in small volumes is due to chiral symmetry restoration, and reproduce this result in our 
calculations. 

Our RG approach improves our understanding of the mechanisms of finite volume effects 
in QCD, but cannot yet give a model independent extrapolation formula to relate finite 
lattice results to the hadronic world. 

In conclusion, we have discussed effects of quark boundary conditions on the finite volume 
shifts of the pion mass. We used the framework of an RG treatment of the quark-meson 
model to offer a possible mechanism which accounts for quark effects. Our approach shows 
the importance of the fermionic boundary conditions for the pion mass and the pion de- 
cay constant. The differences between the results for periodic and anti-periodic boundary 
conditions increase for increasing pion mass and increasing ratio T/L. Our analysis agrees 
qualitatively with the observations from lattice QCD, in regards to the dependence on quark 
boundary conditions as well as in regards to an apparent drop of the pion mass in finite vol- 
ume. We find convergence of our results to those of chiral perturbation theory calculations 
for large pion masses and large volumes, where quark effects are not important. 
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